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ABSTRACT 


The  proposed  greater  reliance  upon  airlifting  military  forces 
demands  that  cargo  loading  time  be  minimized  while  utilization  of 
aircraft  cargo  compartment  space  is  maximized.      Two  loading 
algorithms  have  been  developed  with  these  goals  in  mind  -  a  two 
dimensional  one  for  loading  cargo  where  all  items  must  be  placed 
on  the  floor,    and  a  three-dimensional  one  for  cargo  which  can  be 
stacked.     The  three-dimensional  algorithm  consists  of  the  two- 
dimensional  algorithm  and  a    special  stacking  algorithm.     Tests 
using  randomly  generated  three-dimensional  cargo  lists  indicate 
that  90%  area  efficiencies  for  the  two-dimensional  and  80%  volumetric 
efficiencies  for  the  three-dimensional  algorithm  are  possible.      These 
algorithms  were  designed  for  either  hand  calculations  or  computer 
calculations. 
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I.     INTRODUCTION 

A.      BACKGROUND 

Analyses  of  transportation  and  storage  problems  have  led  to  the 
development  of  many  computer  algorithms  for  the  simulation  of 
loading  cargo  into  containers.        The  usual  objective  of  such  simula- 
tions is  to  determine  the  number  of  containers   required  for  a  given 
list  of  cargo.     This  is  vital  information  in  the  analysis  of  container 
dimensions,    composition  of  transportation  fleets,    etc.  ;  and  computer 
simulation  makes  it  relatively  easy  to  perform  the  necessary 
parametric  studies. 

Loading  simulations  are  particularly  useful  in  military  logistics 
problems.     Accurate  determination  of  the  number  of  ship  and  air- 
craft sorties  required  for  a  given  logistic  operation  allows  a  mean- 
ingful trade-off  to  be  made  between  the  time     required  for  the 
operation  and  the  number  of  transport  vehicles  to  be  assigned. 
Commercial  problems  also  have  this  same  trade-off  situation,    but 
simulation  may  not  be  as  necessary  because  past  experience  with 
similar  loading  situations  often  provides  the  needed  information. 

Military  planners  frequently  need  sortie  data  for  loading  situations 
which  exist  only  on  paper.      Proposed  new  transport  vehicles  and 
transported  vehicles  have  increased  the  use  of  loading  simulations 
in  order  that  new  logistic  situations  may  be  evaluated.      For  this 
reason,    the  majority  of  loading  simulations  have  been  conducted  by 
military  transportation  agencies  and  their  civilian  contractors  - 


The  word  "con  ainer",    when  used  in  this  paper,    refers  to  any- 
thing which  holds  goods,    whether  warehouse,    parking  lot,    or  the 
cargo  carrying  section  of  any  vehicle. 


Boeing,    Douglas,    Lockheed,    RAC,    Rand,    Mitre,    and  others  -  to 

solve  military  problems;  but  the  methods  used  are  also  applicable 

I 
to  ^commercial  problems. 

The  most  sophisticated  loading  algorithms,    such  as  those  com- 
pared by  Eastman  and  Holladay  [l]  attempt  to  fit  the  cargo  into  the 
containers  in  much  the  same  fashion  as  loading  personnel  do.      The 
simulation  results  compare  closely  with  the  results  obtained  by 
loadmasters  in  the  field,    particularly  when  the  cargo  consists  of 
vehicles  and  large  pallets. 

B.   THE  NEED  FOR  BETTER  LOADING  ALGORITHMS 

There  are  two  shortcomings  of  these  sophisticated  algorithms 
which  limit  their  applicability  to  some  future  loading  problems. 
Nearly  all  of  tlu  m  consider  only  the  two-dimensional  problem  of 
loading  the  container  floor  area.      This  was  due  primarily  to  their 
development  for  vehicle  airlift  problems.     The  sizes  and  weights  of 
the  transported  vehicles  and  the  low  heights  of  the  operational  air- 
craft cargo  compartments  made  stacking  infeasible  for  the  problems 
of  major  interest.      Furthermore,    the  algorithms  usually  attempt  to 
predict,    rather  than  improve,    the  performance  of  loadmasters  [1]. 
Most  studies  assume  that  the  loadmaster's   role  in  measuring  and 
fitting  cargo  into  vacant  spaces  is  indispensable  and  that  his  methods 
leave  little  room  for  improvement.     These  assumptions  become 
more  questionable  as  the  size  of  the  airlifted  force  increases. 

The  advent  of  the  jumbo  aircraft,    particularly  the  C-5,    has 
placed  new  emphasis  on  loading  algorithms.     The  proposed  greater 
reliance  on  airlifting  military  forces  in  response  to  threats  has 
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generated  a  need  for  improvement  in  the  algorithms  in  order  to 
make  airlifting  as  effective  as  possible. 

Two  important  goals  of  effective  airlifting  are  rapid  loading  and 
efficient  use  of  aircraft  capacity.     The  present  methods  of  actually 
loading  aircraft  cause  these  goals  to  conflict.     A  loadmaster  usually 
obtains  decreased  container  capacity  utilization  as  the  time  allowed 
for  loading  decreases,    ceteris  paribus.      There  is  an  excellent  chance 
that  computers  with  efficient  loading  algorithms  can  achieve  both 
goals  much  better  than  unaided  human  loadmasters,    particularly 
when  large  items  are  to  be  loaded. 

Decreased  loading  time  could  be  achieved  by  using  loading  al- 
gorithms to  provide  computer  printouts  of  instructions  to  loadmasters 
detailing  exactly  where  each  item  is  to  be  located  in  each  aircraft. 
This  would  all  but  eliminate  the  time-consuming  trial  and  error 
loading  techniques  presently  employed.  Even  when  late  arrival  of 
cargo  or  aircraft  make  the  current  instructions  useless,  a  nearby 
computer  terminal  could  be  used  to  produce  new  instructions   rapidly. 

In  addition  to  speeding  the  loading  times,    the  loading  algorithms 
should  utilize  container  capacities  as  efficiently  as  practicable, 
preferably  surpassing  the  present  trial  and  error  methods. 

Some  possible  benefits  of  successfully  developing  computerized 
loading  instructions  include: 

(1)  The  time  to  airlift  a  given  military  force  would  be  re- 
duced because  of  decreased  loading  time. 

(2)  The  number  of  aircraft  required  for  a  given  airlift 

capability  would  be   reduced,    lowering  total  procure- 
ment and  operating  costs. 
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(3)    The  congestion  at  the  origin,    enroute,    and  at  destination 
airfields  would  be  reduced  because  of  fewer  sorties  re- 
quired for  a  given  airlift  operation.     This  should  further 
reduce  the  time  and  cost  of  an  operation. 

For  effective  competition  with  loadmasters,    three-dimensional 
loading  algorithms  will  have  to  receive  much  greater  emphasis  than 
they  have  in  the  past.      Computerized  instructions  are  needed  for 
stacking  and  loading  not  only  boxes  and  crates  but  also  vehicles. 

The  greater  cargo  compartment  height  of  the   C-5  has  aroused 
interest  in  stacking  the  smaller  and  lighter  military  vehicles,    such 
as  jeeps,    trailers,    mechanical  mules,    etc.  ,    by  designing  them  with 
lower  profiles  and  using  racks  and  loading  frames  which  may  be 
loaded  and  unloaded  easily.     The  spac        n  an  aircraft  above  most  of 
the  loaded  vehicles  is  presently  rarely  used;  and,    as  a  consequence, 
the  binding  constraint  on  container  capacity  is  usually  loaded  floor 
area,    rather  than  volume  or  weight  [l]. 

Attempts  to  produce  computerized  loading  instructions  could  be 
beneficial  even  if  they  are  not  completely  successful.     It  might  be 
feasible  to  computerize  only  the  loading  of  the  larger  items;  even 
this  would  greatly  simplify  the  loadmaster's  task.      Perhaps  only  the 
two-dimensional  problems  are  appropriate  for  computerization, 
which  would  leave  the  stacking  decisions  to  the  loadmaster.      The 
search  for  efficient  algo-ithms  might  result  in  some  new  ground 
rules  which  would  give  loadmasters  a  better  method  of  practicing 
their  art. 
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C.      RELIANCE  ON  HEURISTIC  METHODS 

The  major  problem  in  developing  efficient  loading  algorithms  is 

I 

the  non-applicability  of  existing  mathematical  programming  tech- 
niques for  efficiently  utilizing  two-  and  three-dimensional  space. 
Knapsack  problem  solution  techniques  must  be  allowed  to  choose  how 
many  of  each  size  of  items  are  loaded,    thus  leaving  some  items  on  a 
given  cargo  list  unloaded.      Moreover,    the  best  two-dimensional  knap- 
sack problem  solution  technique  requires  that  the  container  floor  be 
divided  into  rectangles  of  known  dimensions  before  the  technique  can 
be  applied  [2j« 

Cutting  stock  problem  solution  techniques  have  some  applications 
in  container  utilization  efficiency.     Techniques  for  optimal  solutions 
to  one-dimcnsic  :al  problems  of  fitting  a  given  list  of  cargo  into  the 
minimum  number  of  containers  have  been  developed  [3],     Their  use 
is  limited  to  those  cases  where  assumptions  can  be  made  about  how 
the  second  and  third  dimensions  constrain  loading;  e.  g.  ,    stacking  is 
not  possible,    exactly  two  items  can  always  be  loaded  side  by  side, 
etc.      Two-  and  three-dimensional  cutting  stock  solution  techniques 
will  select  the  best  of  many  patterns  submitted  for  consideration, 
but  the  formulation  of  patterns  is  an  intractable  problem  for  even  a 
small  number  of  different  object  sizes  [4]. 

Most  airlifts  have  enough  different  sizes  of  items  to  be  loaded  to 
easily  violate   some  of  the  assumptions  which  must  be  made  before 
present  mathematical  programming  techniques  can  be  used  to  mini- 
mize the  number  of  aircraft  sorties  required.     It  is  possible  that 
new  mathematical  programming  techniques  will  be  developed  to 
maximize  utilization  of  two-  and  three-dimensional  space  with  less 
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restrictive  assumptions,    but  the  likelihood  of  success  in  the  near 
future  seems  to  be  very  low. 

The  most  promising  method  of  developing  efficient  loading 
algorithms  seems  to  be  to  study  the  effects  of  many  collections  of 
loading  decision  rules  in  order  to  ferret  out  those  which  lead  to 
highest  space  utilization.     Attempting  to  achieve  the  highest  possible 
efficiency  could  become  an  endless  task;  a  more  reasonable  goal 
would  be  to  surpass  the  current  efficiency  of  loadmasters. 

The  many  alternatives  an  algorithm  can  take  when  another  item 
is  to  be  loaded  could  eventually  be  separated  into  three  groups:     those 
most  likely  to  increase  efficiency,    those  most  likely  to  decrease 
efficiency,    and  those  whose  effects  can  not  be  safely  predicted.      This 
triage  alone  would  be  a  major  step  forward  in  the  quest  for  efficient 
algorithms. 

D.      A  STANDARD  FOR  LOADING  ALGORITHM  COMPARISONS 

Published  algorithms  with  proven  success  in  predicting  the  load- 
master's  efficiency  were  sought  for  standards  of  comparison.     Only 
one  was  found;  it  was  part  of  a  two-dimensional  loading  model  devel- 
oped by  the  Institute  for  Defense  Analyses  (IDA)  to  determine  the 
number  of  aircraft  sorties  needed  to  haul  a  given  list  of  cargo  [l]. 
The  algoritr         known  as  LOAD  VEHICLES,    loads  rectangles 
(vehicles)  into  a  larger  rectangle  (cargo  compartment  floor).     It  is 
frequently  used  for  airlift  simulations  where  the  "highest  area  effi- 
ciency practicable"  is  desired.     Several  other  loading  models  use 
variations  of  this  algorithm.      Boeing  Aircraft's  SLAM  program, 
whose  efficiency  is  compared  with  the  IDA  model  in  Ref.    1,    is  an 
example. 
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H.       OBJECTIVES  AND  SCOPE 

The  objectives  of  this  study  were  to  develop  a  volumetrically 
efficient  three-dimensional  loading  algorithm  and  an  area  efficient 
two-dimensional  loading  algorithm.     Such  algorithms  are  the  first 
step  in  computerizing  loading  instructions.      It  was  considered  desir- 
able,   but  not  essential,    that  any  algorithms  developed  would  utilize 
container  space  as  well  as  or  better  than  the  average  loadmaster. 

It  was  also  considered  desirable,    but  not  necessary,    that  the 
algorithms  developed  be  simple  enough  for  use  without  a  computer 
and  that  integer  and  linear  programming  be  avoided  if  no  significant 
efficiency  would  be  lost  thereby.     This  was  accomplished,    greatly 
reducing  computer  time,    and  it  made  the  algorithms  usable  for  hand 
calculations  by  loading  per*  onnel. 

The  three-dimensional  problem  was  conceived  as  the  filling  of 
rectangular  solid  containers  with  rectangular  solid  objects  so  as  to 
minimize  the  number  of  containers  required  to  hold  any  given  number 
of  objects  of  many  different  sizes.     Similarly,    the  two-dimensional 
problem  involved  only  the  rectangular  floors  of  the  containers  and 
the  rectangular  bases  of  the  objects  to  be  loaded. 

Container  door  dimensions  were  not  considered.      For  simplicity, 
each  piece  of  cargo  was  assumed  to  be  marked  "THIS  END  UP". 

Cargo  weight  was  not  considered  in  the  algorithms  tested.     The 
containers  were  assumed  to  be  stror      enough  to  support  any  weight 
placed  anywhere.      Center  of  gravity  movement  was  disregarded. 
Each  cargo  item  was  assumed  to  be  strong  enough  to  support  what- 
ever other  items  that  might  be  placed  or  top  of  it. 
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The  algorithms,    as  presented,    do  not  print  out  the  locations  of 

the  cargo  items;  but  simple  modifications  to  the  computer  program 

i 

on  pages  40  through  43   allow  this. 
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III.       THE  SYNTHETIC  TEST  LOADS 

The  desire  to  efficiently  load  all  the  items  from  a  list  containing 
different  quantities  of  each  of  many  different  sizes  of  items  required 
some  method  of  selecting  cargo  dimensions  and  container  dimensions. 
A  group  of  dissimilar  cargo  lists  was  selected  for  testing  proposed 
algorithms. 

The  unit  of  length  measurement  was  set  at  six  inches,    and  only 
integer  values  were  used  throughout  the  study.     This  was  to  restrict 
the  number  of  different  sets  of  dimensions  to  a  more  manageable 
group.     A  larger  length  unit  would  have  reduced  the  group  further, 
but  it  would  also  have  reduced  the  accuracy  with  which  the  cargo 
items  could  be  measured,    since  each  item's  dimensions  are  rounded 
up  to  the  next  integer  unit. 

A  notio'-al  aircraft  cargo  compartment  was  selected  with  a  length 
of  200  units,    width  of  30  units,    and  a  height  of  20  units.     These  dimen- 
sions were  not  changed  throughout  the  study  because  the  desired 
variety  in  the  loads  was  achieved  by  varying  the  list  of  cargo. 

Cargo  lists  were  made  by  computer  generation  of  one  item  at  a 
time  until  a  selected  volume  was  exceeded.     Two  random  selections 
were  made  from  the  integers  1,2,...,  20,    where  each  integer  had  an 
equal  probability  of  being  selected.     The  larger  integer  selected  was 
the  item's  length;  the  smaller,    the  width.     The  height  was   similarly 
selected  from  1,2,...,  10.     This  produced  2,  100  different  sets  of 
dimensions  for  items  which  could  be  on  the  cargo  lists.     The  expected 
values  of  the  length,    width,    and  height  were  13.825,    7.  175,    and  5.5 
units,    respectively. 
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Tables  I  and  II  show  two  distributions  of  generated  cargo  item 
volumes  for  four  approximate  total  cargo  volumes.     Each  table  was 
compiled  by  generating  enough  items  to  exceed  80,  000  cubic  units  of 
volume  and  the  distribution  of  item  volumes  was  recorded.      Next, 
more  items  were  generated  and  their  volumes  added  to  the  80,  000 
until  120,  000  cubic  units  was  exceeded.      The  recording  and  generation 
of  additional  items  were  continued  using  total  volume  increments  of 
40,  000  cubic  units  until  data  was  obtained  for  200,  000  cubic  units. 
These  tables  will  be  used  to  illustrate  the  results  obtained  by  the 
algorithms  presented  below. 

Many  other  cargo  lists  were  generated  and  loaded,    but  the 
algorithms  produced  no  interesting  changes  in  resul'. s  for  the  same 
approximate  total  volumes  of  cargo.     A  few  tests  with  larger  average 
dimensions  for  the  cargo  items  confirmed  intuition  that  loading 
efficiencies  decreased  for  the  two  algorithms  presented  below. 


1  APPROXIMATE 
! TOTAL 
VOLUME 

80,000 

120,000 

160,000 

200,000 

GENERATED 

TOTAL 

VOLUME 

80,629 

120,172 

160,623 

200,487 

TOTAL 

NUMBER 
OF  ITEMS 

133 

206 

285 

361 

NUMBER  OF 
ITEMS  IN  EACH 
VOLUME  RANGE 

CO 
H 
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CJ 
M 
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M 
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0  -  49 

11 

16 

24 

37 

50  -  sy 

13 

20 

28 

31 

100  -  149 

11 

17 

23 

27 

150  -  1° 

8 

15 

22 

29 

200  -  299 

20 

28 

35 

43 

300  -  399 

8 

15 

24 

30 

400  -  599 

15 

26 

32 

41 

600  -  799 

12 

17 

25 

36 

800  -  999 

12 

16 

23 

24 

1,000  -  1,499 

8 

15 

25 

34 

1,500  -  1,999 

5 

9 

11 

16 

2,000  -  4,000 

10 

12 

13 

13 

Table  I.    Sample  distribution  of  cargo  items  by  volume 
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APPROXIMATE 

TOTAL 

VOLUME 

80,000 

120,000 

160,000 

200,000 

GENERATED 

TOTAL 

VOLUME 

81,589 

120,103 

160,138 

200,005 

TOTAL 

NUMBER 
OF  ITEMS 

162 

220 

292 

374 

NUMBER  OF 
ITEMS  IN  EACH 
VOLUME  RANGE 

to 

H 
M 

O 

M 
CQ 
!=> 

& 
M 

o 
> 

0-49 

20 

25 

30 

46 

50  -  99 

16 

24 

39 

47 

100  -  ."  \9 

17 

24 

30 

37 

150  -  199 

18 

22 

26 

28 

200  -  299 

15 

18 

22 

29 

300  -  399 

12 

14 

21 

29 

400  -  599 

18 

27 

34 

43 

600  -  799 

16 

20 

23 

32 

800  -  999 

5 

9 

14 

19 

1,000  -  1,499 

12 

16 

23 

28 

1,500  -  1,999 

6 

8 

16 

20 

2,000  -  4,000 

7 

13 

14 

16 

Table  II.    Sample  distribution  of  cargo  items  by  volume 
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IV.       THE  STACKING  ALGORITHM 

A.      GENERAL 

The  stacking  algorithm  to  be  presented  below  was  selected  as 
the  simplest  method  of  achieving  a  good  three-dimensional  loading 
algorithm.      It  reduces  the  three-dimensional  problem  to  a  two- 
dimensional  problem  by  loading  the  cargo  onto  notional  pallets  which 
must  then  be  loaded  into  the  containers  by  any  floor  area  loading 
algorithm. 

The  base  of  each  notional  pallet  is  the  largest  cargo  item  on  that 
particular  pallet,    rather  than  literally  a  metal  or  wooden  platform 
under  the  cargo.     Therefore,    each  pallet  base  takes  the  dimensions 
of  the  base  of  its  largest  cargo  item. 

The  objective  of  the  stacking  algorithm  is  to  maximize  stacking 
efficiency  where: 

(1)  STACKING  EFFICIENCY    =    TOTAL  CARGO  VOLUME 


TOTAL  STACKING 

PALI.  FT  HEIGHT 

AREA 

Stacking  height  is  a  constant  equal  to  container  height  throughout  this 

paper.     Thus,    the  algorithm  achieves  its  objective  by  minimizing  the 

total  area  of  the  notional  pallets  it  must  use  for  a  given  volume  of 

cargo  to  be  loaded. 

Volumetric  loading  efficiency  is  the  stacking  efficiency  of  the 

stacking  algorithm  multiplied  by  the  floor  area  loading  efficiency 

obtained  by  the  two-dimensional  algorithm  which  loads  the  notional 

pallets  into  the  aircraft.     The  use  of  a  stacking  algorithm  to  paFctize 
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the  cargo  before  anything  is  placed  on  the  container  floor  simplifies 
efforts  to  improve  volumetric  efficiency  because  such  efforts  can  be 
divided  between  two  paths  which  are  considerably  less  complicated 
than  trying  to  search  for  improvements  in  some  three-dimensional 
algorithm. 

B.      METHOD 

The  first  step  in  the  algorithm  is  to  order  all  of  the  carcio  items 
by  base  size  (area  or  perimeter).     Items  with  the  same  base  size 
are  ordered  by  height.     The  stacking  begins  with  the  largest  item 
being  designated  "Stack  #1";  its  length,   width,    height,    and  base  area 
are  recorded.     The  second  largest  item  is  then  compared  with  the 
top  of  Stack  #1.     The  item  is  stacked  on  #1  if  it  does  not  overhang 
any  side  of  the  top  and  if  it  does  not  cause  the  stack's  height  to  exceed 
the  stacking  height;  otherwise,    it  becomes  "Stack  #2". 

Whenever  an  item  is  loaded  onto  a  stack,    its  base  dimensions 
become  the  new  dimensions  for  the  top  of  the  stack,    and  its  height  is 
subtracted     from  the  stack's  ceiling  clearance  to  obtain  the  new  clear- 
ance.     The  unused  area  on  the  previous  top  of  the  stack  is  used  for  a 
substack.      The  base  of  the  substack  is  the  larger  area  rectangle 
ABGH  or  AJED  of  Figure  1.     Substacking  is  performed  using  the 
same  rules  of  fit  as  stacking.      Items  are  substacked  from  the  cargo 
list  until  either  the  substack  reaches  the  ceiling  or  the  entire  cargo 
list  has  been  scanned  for  items  yet  to  be  loaded  which  will  fit  onto 
the  substack.     Substacks  are  not  numbered.     The  algorithm  "forgets" 
them  after  every  effort  to  fill  them  is  completed. 
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Figure  1.       Sub  stack  bases 


Stacking  then  resumes  with  the  largest  item  2  emaining  on  the 
cargo  list.     An  attempt  is  made  to  stack  the  item  on  the  unfilled  stacks 
in  order  of  stack  serial  numbers,    i.  e.  ,    in  the  order  in  which  the 
stacks  were  created.      If  the  item  fits  one  of  the   stacks,    a  substack  is 
placed  on  that  stack,    if  possible.      If  the  item  can  not  be  placed  on  any 
unfilled  stack,    it  becomes  the  base  for  a  new  stack. 

The  stacking  algorithm  continues  until  every  item  on  the  cargo 
list  is  positioned  in  a  stack  as  the  base,    a  member  of  the  stack  proper, 
or  a  member  of  one  of  the  several  substacks  of  the   stack.      Figure    2 
is  a  schematic  diagram  of  the  stacking  process. 
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Figure  2.    Schematic  diagram  of  the  stacking  process 
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C.      TEST  RESULTS 

The  stacking  algorithm  proved  to  be  quite  efficient.      Table  III 
shows  the  stacking  efficiencies  obtained  for  the  cargo  volumes 
listed  in  Tables  I  and  II.      These  results  include  both  ordering  the 
cargo  items  by  base  area  and  base  perimeter. 


APPROXIMATE 
CARGO  VOLUME 

80,  000 

120, 000 

160,  000 

200,  000 

CARGO  FROM 
TABLE  I 

STACKING 
EFFICIENCY 
(AREA  ORDERING) 

.  908 

.922 

.938 

.  943 

STACKING 
EFFICIENCY 
(PERIM.    ORDERING) 

.  906 

.  927 

.934 

.  941 

CARGO  FROM 
TABLE  II 

STACKING 
EFFICIENCY 
(AREA  ORDERING) 

.  913 

.935 

.945 

.  950 

STACKING 
EFFICIENCY 
(PERIM.    ORDERING) 

.  907 

.918 

.943 

.  950 

Table  III.     Efficiencies  of  the   slacking  algorithm 

The  most  interesting  test  result  shown  in  Table  III  is  the  general 
tendency  of  stacking  efficiency  to  increase  with  increased  volume  of 
cargo.      This  higher  efficiency  with  higher  volume  was  expected 
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because  a  larger  number  of  items   selected  from  a  fixed  number  of 
different  dimensions  should  cause  the  base  area  of  each  stacked  item 
to  more  nearly  cover  the  base  area  of  the  item  below  it  in  the  stack. 

Tests  with  cargo  volumes  between  the  tabulated  volumes  of 
Table  III  showed  that  this   stacking  efficiency  increase  is  not  monotonic 
as  the  table  might  imply.     Additional  cargo  items  were  generated  and 
added  to  those  of  Table  I  to  produce  a  total  cargo  volume  of  600,  000 
cubic  units;  this  volume  was   stacked  with  98.  0%  efficiency  after  area 
ordering. 

This  same  increased  efficiency  phenomenon  obviously  made  the 
use  of  substacking  less  important  as  cargo  volume  increased.     Sub- 
stacking  increased  the  efficiency  by  about  2.  5%  for  approximate  total 
cargo  volumes  of  100,  000  cubic  units,    but  the  increase  was  less  than 
1.  0%  when  total  cargo  volume  exceeded  400,  000  cubic  units.      Sub- 
stacking  thus  appears  to  provide  very  little  extra  efficiency  for  the 
synthetic  test  loads.      It  was  not  deleted  from  the  algorithm  because 
it  can  provide  much  larger  efficiency  increases  in  situations  where 
there  is  a  large  difference  in  the  base  area  of  an  item  and  the  next 
smaller  item  on  the  cargo  list. 

Simulations  were  conducted  with  two  realistic  vehicle  lists  to 
determine  how  much  use  the  stacking  algorithm  could  make  of  the 
space  in  a  C-5  if  stacking  were  feasible.      The  first  list  was  that  for 
the  82nd  Airborne  Division,    published  in  Ref.    1,    containing  1,  573 
vehicles  of  43  different  types.      The   second  list  was  for  an  infantry 
division  and  contained  6,  811   vehicles  of  132  different  types.     The 
first  list  was   stacked  v/ith  80.  0%   efficiency  and  produced  763   stacks. 
When  loaded  by  the  IDA  algorithm,    these   stack,    required  only  38 
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aircraft  as  opposed  to  56  for  the  unstacked  case.     The  second  list 
resulted  in  2,  664  stacks  with  84.  0%  efficiency.      The  IDA  algorithm 
loaded  these  stacks  into  186  aircraft  vice  326  for  the  unstacked  case. 

At  the  end  of  all  tests,  the  question  of  ordering  the  cargo  items 
by  base  area  or  base  perimeter  was  resolved  in  favor  of  base  area. 
Area  ordering  usually,  but  not  always,  produced  a  slightly  higher 
stacking  efficiency  than  that  obtained  by  perimeter  ordering.  Table 
III  is  typical  of  the  results  obtained  for  all  cargo  lists  tested.  Note 
that  the  one  case  where  area  ordering  was  not  superior  is  in  the 
120,  000  cubic  units  column. 

D.      SUMMARY 

The  stacking  algorithm  presented  combines  high  stacking  efficiency 
with  computational  ease.      In  the  tests  it  provided  increased  efficiency 
as  cargo  volume  was  increased.      It  can  be  used  with  any  ;v/o-dimension- 
al  lo:   Hng  algorithm  to  produce  a  three    dimensional  algorithm. 

A  major  obstacle  to  the  use  of  stacking  algorithms  is  the  nature 
of  the  cargo.     They  will  have  limited  applications  for  loading  vehicles 
until  engineering  permits  the  stacking  of  some  vehicles  in  airlift 
situations. 
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V.       THE  LENGTH-MODULAR  ALGORITHM 

A.  I    GENERAL 

The  length -modular  algorithm  is   so  named  because  it  divides 
each  container  into  modules  whose  length     are  that  of  the  longest 
cargo  item  in  each  module  and  whose  widths  are  equal  to  container 
width.      It  is  a  two-dimensional  loading  algorithm  which  attempts  to 
maximize  area  efficiency  where: 

(2)     AREA  EFFICIENCY    =    TOTAL  PALLET  AREA  LOADED 


NUMBER  OF  CONTAINER 

CONTAINERS       X       BASE  AREA 
REQUIRED 

As  in  equation  (1),  the  pallet  of  equation  (2)  is  notional.  There  may 
or  may  not  liberally  be  a  wooden  or  metal  platform  under  everything 
which  covers  some  part  of  the  container  floor. 

B.      METHOD 

The  algorithm  loads  one  container  at  a  time.      It  first  finds  the 
length  of  the  longest  cargo  item  and  uses  that  as  the  length  of  the 
first  module.     Whenever  a  module  is   created,    it  is  partitioned  into 
three  rectangles  as   shown  in  Figure  3.     Rectangle  A  is  always  com- 
pletely covered  with  cargo,    its  length  is  always  equal  to  that  of  the 
module,    and  its  wdth  is  initially  zero.     Rectangle  B,    known  as  a 
submodule,    is  always  empty  and  it  initially  covers  the  entire  module. 
Rectangle   C,    known  as  a  lateral,    maybe  empty,    partially  filled,    or 
completely  filled;  its  initial  dimensions  are  both  zero.      The  lengths 
of  rectangles  A  and  B  are  measured  in  the  same  direction  as  the 

module  length,   but  the  length  of  rectangle  C  is  always  its  longer 
dimension. 
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MODULE  AND  CONTAINER  WIDTH 


Figure  3.       Partitions  of  a  module 


The  first  step  in  filling  the  module  is  to  load  the  widest  item  with 
the  same  length  as  the  module.     The  widths  of  rectangles  A  and  B  are 
respectively  increased  and  decreased  by  the  item's  width.     The  al- 
gorithm continues  loading  the  widest  items  of  that  length  until  the 
width  of  rectangi:    B  is  zero,    all  items  of  that  length  are  loaded,    or 
all  items  yet  to  be  loaded  of  that  length  are  wider  than  the  width  of 
rectangle  B.     Whenever  the  first  case  occurs,    rectangle  A  occupies 
the  entire  module,    the  algorithm  is  finished  with  that  module,    and  a 
new  module  is  started. 

If  the  first  case  does  not  occur,    the  algorithm  finds  the  longest 
possible  item  which  will  fit  into  rectangle  B.     If  none  will  fit,    the 
algorithm  is  finished  with  that  module  and  starts  a  new  one.      If  some 
item  does  fit  into  rectangle  B,    say  vehicle  X,    the   rectangle's  length 
is  decreased  to  the  length  of  vehicle  X.     The  area  lost  by  rectangle  B 
because  of  this  decrease  becomes  the  area  of  rectangle   C.      Note  that 
vehicle  X  has  not  been  loaded  yet.     Rectangle  C,    the  lateral,    is  then 
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loaded  with  the  longest  item  possible,    then  the  next  longest  possible, 
etc.  ,    until  no  items  on  the  cargo  list  will  fit  lengthwise  into  the  un- 
loaded length  of  the  rectangle.     Rectangle   C  is  loaded  before  rectangle 
B  in  order  that  the  same  computer  routine  which  loads  an  empty- 
module  can  load  an  empty  submodule  without  having  to  return  to 
rectangle   C  or  remember  its  dimensions. 

Figure  4  shows  a  sample  filled  lateral.      Note  that  item  3  could 
have  been  rotated  90  degrees  to  make  the  unused  length  of  the  lateral 
longer.     This   rotation  was  found  to  be  unproductive,    as  a  rule,    be- 
cause laterals  are  almost  always  extremely  narrow,    on  the  order  of 
two  or  three  units  wide  for  the  synthetic  loads.      They  would  also 
normally  be  too  narrow  to  hold  the  shortest  vehicle  in  an  actual  load. 


Figure  4.       Sample  load  in  a  lateral 

The  algorithm  then  prepares  to  load  rectangle  B  with  the  longest 
item  possible.      If  this  item  is  not  vehicle  X,    which  might  have  been 
loaded  into  rectangle   C,    or  one  of  the  same  length,    the  length  of 
rectangle  B  is  further  decreased  to  the  length  of  this  new  longest  item. 
The  area  lost  in  rectangle  B  is  again  added  to  rectangle  C,    but  no 
effort  is  made  to  fill  this  area  because  this  situation  rarely  occurs 
during  a  load  of  many  items. 
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Rectangle  B,    the  submodule,    is  then  loaded  in  the  same  manner 

as  the  original  module,    i.e.  ,    it  is  partitioned  into  three  rectangles, 

i 
A',    B',    C',    which  are  loaded  in  the  same  manner  as  the  original  A, 

B,    C.      This  partitioning  of  submodule s  continues  until  some  sub- 
module  of  a  submodule  is  too  small  for  any  unloaded  item  on  the 
cargo  list.      The  maximum  number  of  partitionings  of  a  module  and 
its  submodules  is  the  module's  length  in  integer  units;  the  minimum 
is  one,    regardless  of  length.      Four  or  more  partitionings  were  ex- 
tremely rare  for  the  synthetic  test  loads. 

Figure  5  illustrates  how  a  module  might  look  when  the  algorithm 
is  finished  with  it.      The  rectangles  A,   A',   A",   A'"  are  completely 
filled.     The  rectangles   C,    C,    C"  may  each  be  filled,    partially  filled, 
or  empty.     Rectangle  B'"  is  empty,    but  its  area  could  be  zero. 


A 

A1 

A" 

A"' 

B1" 

C" 

C 

C 

Figure  5.       Sample  final  partitioning  of  a  module 

The  second  and  su-         ij.ng  modules  in  a  container  are  given  a 
length  equal  to  the  longest  cargo  item  remaining  on  the  list  which  is 
not  longer  than  the  remaining  length  of  the  container.     This  makes 
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each  module  no  longer  than  any  created  before  it  in  that  container. 
The  process  of  creating  and  loading  modules  continues  until  some 
module  is  created  which  is  too  short  for  all  of  the  remaining  items 
on  the  cargo  list.      If  all  items  are  loaded,    the  algorithm  is  finished; 
otherwise,    it  starts  with  a  new  container.      Figure  6  is  a  schematic 
diagram  of  the  length-modular  algorithm. 

C.      TEST  RESULTS 

The  length-modular  algorithm  was  tested  in  conjunction  with  the 
stacking  algorithm  for  many  synthetic  test  loads.     The  stacking  al- 
gorithm was  usually  applied  first,    and  then  the  length-modular 
algorithm  loaded  the  stacks  into  the  containers. 

Table  IV  shows  the  results  of  applying  the  length- modular  algo- 
rithm to  the  pallets  which  were  stacked  and  then  tabulat  .  1  in  Table  III. 
The  results  of  loading  the  same  pallets  with  the  IDA  algorithm  are 
a^o  presented  in  Table  IV  for  comparison.     Area  efficiency  was 
computed  in  each  case  by  treating  the  last  loaded  container's  length 
equal  only  to  the  loaded  length  of  the  container.     This  permitted  a 
more  meaningful  comparison  of  the  efficiencies  to  be  made,    since  no 
more  than  three  containers  were  ever  required  for  the  volumes  in 
the  preceding  tables. 

Area  efficiency  for  the  length-modular  algorithm  was  found  to 
increase  with  increased  cargo  volume  in  a  manner  similar  to  that  of 
the  stacking  algorithm.     A  larger  number  of  items   selected  from  the 
210  different  sets  of  base  dimensions  usually  enabled  the  algorithm 
to  find  a  better  fitting  cargo  item  for  a  particular  vacant  space  than 
when  the  number  of  items  to  be  loaded  was  fewer. 
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Figure  6.    Schematic  diagram  of  the  length-modular  algorithm 
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APPROXIMATE 
TOTAL  CARGO 
VOLUME 

80, 000 

120,  000 

160,  000 

200,  000 

TABLE  I  CARGO 

STACKING 
EFFICIENCY 
(AREA  ORDERING) 

.  908 

.  922 

.  938 

.  943 

AREA  EFFICIENCY 
WHEN  THESE  STACKS 
WERE  LOADED  BY 
THE  L-M  ALGORITHM 

.  912 

.  935 

.  952 

.965 

AREA  EFFICIENCY 
WHEN  SAME  STACKS 
WERE  LOADED  BY 
THE  IDA  ALGORITHM 

.790 

.  780 

.766 

.  767 

TABLE  II  CARGO 

STACKING 
EFFICIENCY 
(AREA  ORDERING 

.913 

.  935 

.  945 

.  950 

AREA  EFFICIENCY 
WHEN  THESE  STACKS 
WERE  LOADED  BY 
THE  L-M  ALGORITHM 

.903 

.  952 

.  943 

.960 

AREA  EFFICIENCY 
WHEN  SAME  STACKS 
WERE  LOADED  BY 
THE  IDA  ALGORITHM 

.  767 

.  775 

.  781 

.  774 

Table  IV:      Area  efficiency  comparison  of  the  length- modular  and 
IDA  algorithms 
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Three  other  permutations  of  modular  loading  were  tested.     The 
first  might  be  called  "width- modular"  because  module  length  was  set 
equal  to  the  width  instead  of  the  length  of  the  largest  remaining  item. 
The  other  two  permutations  were  length-modular  and  width- modular 
with  module  width  set  equal  to  the  container's  length  instead  of  width. 
These  permutations  gave  much  poorer  results  than  the  original 
method. 

The  same  two  vehicle  lists  discussed  in  Section  IV  were  loaded 
unstacked  by  the  length-modular  algorithm,    and  comparisons  were 
made  with  the  IDA  algorithm.     The  smaller  list  was  loaded  into     54 
C-5's  vice  56  for  the  IDA  model.     The  larger  list  required  only  310 
aircraft  vice  the  IDA  model's  326. 

It  was  interesting  to  note  that  the  IDA  model  gave  area  efficiencies 
within  +  3.  2%  of  78.  0%  throughout  about  40  synthetic  loads  of  ten 
different  volumes.     However,    when  it  loaded  the  two  lists  of  vehicles, 
area  efficiency  increased  to  about  84%  for  both  lists.     Its  performance 
did  not  ever  equal  that  of  the  length-modular  algorithm  for  any  of  the 
tests  conducted,    but  there  could  be  cases  where  it  would  be  superior. 

A  few  tests  were  conducted  where  the  length- modular  algorithm 
loaded  the  container  floor  first,    and  then  the  stacking  algirithm  loaded 
vertically  upon  those  items  which  covered  the  floor.      The  area  efficien- 
cy thus  obtained  was  very  high,    never  less  than  96.  7%,    but  the  stacking 
efficiency  was  so  degraded  that  volumetric  efficiency  was  always 
15-    0%  lower  than  when  the  same  load  was  stacked  before  the  container 
floor  was  covered. 
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D.      SUMMARY 

The  length-modular  algorithm  has  been  shown  to  provide  excellent 
area  efficiency  for  some  loading  situations.     It  is  simple  in  method, 
although  the  computer  program  on  pages  40  through  43   is  somewhat 
complicated  by  steps  to  simplify  record  keeping  and  reduce  computer 
time.     There  are  several  places  where  the  algorithm  as  listed 
sacrifices  area  efficiency  in  order  to  save  time.      It  generally  provides 
increased  efficiency  with  increased  total  volume  of  cargo. 

It  should  be  noted  that  the  modules  in  any  loaded  container  can  be 
repositioned  to  move  the  container's  center  of  gravity  longitudinally, 
and  items  within  each  module  can  be  moved  in  several  ways  to  move 
the  center  of  gravity  laterally.      This  feature  of  modular  loading 
facilitates  the  algorithm's  proposed  use  for  coi  lputerizing  loading 
plans  for  aircraft. 
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VI.   AREAS  FOR  FURTHER  STUDY 

The  particular  rules  and  methods  used  in  this  study  are  only  a 
minute  part  of  what  could  be  considered.      The  high  efficiencies  of 
the  algorithms  presented  will  not  permit  major  increases,    but  some 
worthwhile  increases  in  efficiency  might  be  easily  discovered  in  both 
algorithms.      The  areas  suggested  below  for  further  study  are  only  a 
few  which  might  be  promising  for  increased  efficiency  and  inclusion 
of  aspects  such  as  weight  and  center  of  gravity. 

A.  METHODS  WITH  "THIS  END  UP"  ASSUMPTION  REMOVED 
Since  many  cargo  items  may  be  loaded  with  any  of  its  three  axes 

vertical,    it  would  be  useful  to  know  which  axis  should  be  placed 
vertically  when  the  stacking  algorithm  is  given  a  choice.     Simple 
rules,    such  as  prescribing  the  longest,    shortest,    or  middle  length 
axis,    might  be  found  to  yield  the  highest  volumetric  efficiency.      More 
complex  rules,    which  select  a  different  axis  for  differed  1  stack 
clearances  or  other  stacking  parameters,    might  be  necessary. 

B.  PRE-STACKING 

There  are  countless  ways  in  which  many  items  of  one  or  more 
common  dimensions  might  be  combined  into  a  rectangular  solid 
having  little  or  no  wasted  space.     Such  a  solid  would  then  be  stacked 
as  one  item.      This  type  of  cargo  list  consolidation  befo;  .    the  stacking 
algorithm  is  applied  might  have  surprising  advantages  in  efficiency 
and  speed. 
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C.  HEIGHT -MODULAR  STACKING 

It  would  be  interesting  to  know  what  might  be  done  by  having  a 
two-dimensional  loading  algorithm  load  modules  with  only  items  of 
particular  height  ranges  and  then  stack  the  modules  in  various  ways. 
Such  modules  would  not  necessarily  cover  the  container  floor. 

D.  WEIGHT  CONSIDERATIONS 

Modifications  to  the  two  algorithms  presented  could  allow  them 
to  consider  center  of  gravity  constraints  and  total  weight  applied  to 
any  part  of  the  container  floor.     This  will  be  necessary  before  com- 
puterized loading  instructions  can  become  a  reality  for  aircraft.     The 
modifications  might  not  significantly  decrease  the  efficiencies  of  the 
algorithms. 

E.  NON-RECTANGULAR  CONTAINERS 

Methods  for  loading  non- rectangular  containers  with  rectangular 
cargo  have  received  even  less  attention  than  the  rectangular  container 
case.     Aerodynamic,    hydrodynamic,    and  other  engineering  considera- 
tions dictate  that  many  containers  take  on  shapes  which  will  always 
result  in  some  wasted  space  for  any  realistic  non-fluid  load. 
Minimization  of  the  wasted  space  would  be  a  real  challenge. 
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VII.       CONCLUSIONS  AND  RECOMMENDATIONS 

The  two  algorithm?;  preserted  have  demonstrated  high  volumetric 
and  area  efficiencies  for  loading  a  large  number  of  items  of  many- 
different  sizes.     They  are  an  important  first  step  in  computerizing 
the  loading  of  aircraft  and  other  containers.     The  length-modular 
algorithm  is  the  more  important  of  the  two  for  major  airlifting 
problems  because  they  presently  involve  cargo  which  permits  little 
stacking. 

Much  remains  to  be  done  before  computerized  loading  can  become 
a  reality.      Loading  algorithms  must  be  able  to  consider  each  item's 
weight  and  fragility,    as  well  as  its  effect  upon  the  container's  center 
of  gravity.      Methods  for  utilizing  container  space  more  efficiently 
should  be  sought,    but  the  algorithms  presented  here   should  be  good 
enough  for  come  of  the  interim  work  necessary  for  development  of 
computerized  loading  systems. 

At  this  stage,    a  logical  next  step  would  be  to  see  how  well  the 
length-modular  algorithm  can  compete  with  loadmasters  in  two- 
dimensional  loading  of  vehicles  without  weight  constraints.      If  the 
IDA  loading  algorithm  is  truly  an  accurate  predictor  of  human  loading 
efficiency,    then  chances  are  excellent  that  the  length-modular  al- 
gorithm can  make  computerized  aircraft  loading  a  reality. 
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COMPUTER  PROGRAM 


C  THE  NPS  LOAD  MODEL 

C  DETERMINES  THE  NUMBER  OF  CONTAINERS  TO  HOLD  THE  INPUT 

C  LIST  OF  ITEMS 

C  LIMITATIONS:  7,000  ITEMS  OF  LENGTH  NOT  GREATER  THAN  100 
C  UNITS,  WIDTH  OF  40  UNITS,  HEIGHT  OF  99  UNITS, 

C  UP  TO  300  DIFFERENT  TYPES  OF  ITEMS  ALLOWED 

C  UP  TO  999  CONTAINERS  ALLOWED,  ALL  MUST  HAVE 

C  SAME  DIMENSIONS,  LENGTH  NOT  GREATER 

C  THAN  999  UNITS,  WIDTH  99,  HEIGHT  99 

C  FIRST  DATA  CARD  HAS  NUMBER  OF  ITEM  TYPES  AND  CONTAINER 

C  LENGTH,  WIDTH,  AND  HEIGHT  IN  FOUR  TEN  COLUMN  FIELDS 

C  ALL  OTHER  DATA  CARDS  HAVE  NUMBER  OF  ITEMS  OF  THAT  TYPE 

C  AND  THE  LENGTH,  WIDTH,  AND  HEIGHT  OF  THAT  TV  E  IN  FOUR 

C  TEN  COLUMN  FIELDS 


;   INITIALIZING 

DIMENSION  I!  (300),  IWOOO),  IHC300),  ML(1600). 
lMW(lf'"»0),  MH(1600),  N(300),  LL(I600J,  LW(1600), 
2KU1600),  KWC1600J,  IS(300),  NE(1600),  LH(1600), 
3MTX( 100,41  ) 

DO  600  K  =  1,  100 

DO  600  L  =  1,  41 
600   MTX(K,L)  =  0 

NV  =  0 

JP  =  0 

IB  =  0 

READ(5,91)  MI,  IAL,  IAW,  IAH 
91   FORMAT(4I10) 

DO  954  K  =  1,  Nil 

READ   (5,91)  M,  I RL ,  I RW ,  IRH 

NV  =  NV  +  M  *  IRL  *  IRW  *  IRH 

IB  =  IB  +  M 

JP  =  JP  +  1 

IL(JP)  =  IRL 

IW(JP)  =  IRW 

IH(JP)  =  IRH 

N(JP)  =  M 
954   IS( JP)  =  JP 

:   COMMENCE  SORTING 
JD  =  JP  -  1 
DO  100  J  =  1,  JD 
M  =  IS(J) 
IQ  =  J 
JW  =  J  +  1 
DO  200  I  =  JW,  JP 
MG  =  IS( I ) 
IF  (  IL(M)  *  IW(M)  -  IL(MG)  *  IW(MG)  )  40,  50,  200 


40 


50   IF  (  IH(M)  .GE.  IH(MG)  )  GO  TO  200 
40   M  =  MG 

IQ  =  I 
200   CONTINUE 

IDUM  =  N(IQ) 

N(IO)  =  N(J) 

NU)  =  IDUM 

IS( 10)  =  IS( J) 

IS(J)  =  M 
100   CONTINUE 

;   COMMFNCE  STACKING 

DO  10  K  =  1,  IB 

ML(K)  =  100 

MW(K)  =  40 

LH(K)  =  IAH 
10   MH(K)  =  IAH 

LA  =  0 

M  =  0 

L  =  1 

DO  400  MG  =  1,  JP 

IF  (  N(MG)  .EQ.  0  )  GO  TO  400 

IQ  =  N(MG) 

DO  400  J  =  1,  IQ 

JE  =  M  +  1 

I  =  IS(MG) 

DO  300  JZ  =  L,  JE 

K  =  JZ 

IF  (  IL (  I )  .LE.  ML(K)  .AND.  IW(I)  .LE.  MW(K)  .AND. 
3IHU  )  .LE.  MH(K)  )  GO  TO  56 
300   CONT  MUE 
56   IF  (  MH(K)  .EQ.  IAH)  GO  TO  60 

LH(K)  =  MH(K) 

IF  (  (  MW(K)  *  (  ML(K)  -  IL(I)  )  )  .GE.  (  ML ( K )  * 
4(  MW(K)  -  IV'(  I  )  )  )  )  GO  TO  116 

LL(K)  =  ML(K) 

LW(K)  =  MW(K)  -  IW(I) 

GO  TH  121 
116   LL(K)  =  MW(K) 

LW(K)  =  ML(K)  -  IL ( I ) 

IF  (  LL(K)  .GE.  LW(K)  )  GO  TO  121 

IDUM  =  LL(K) 

LL(K)  =  LW(K) 

LW(K)  =  IDUM 
121   DO  440  IC  =  MG,  JP 
128   IF  (  N(IC)  .EQ.  0  )  GO  TO  440 

IF  (  LH(K)  .EQ.  0  )  GO  TO  60 

IE  =  IS( IC) 

IF  (  IL(IE)  .GT.  LL(K)  .OR.  IW(IE)  .GT.  LW(K)  .OR 
5IHCIE)  .GT.  LH(K)  )  GO  TO  440 

N(IC)  =  N( IC)  -  1 

LL(K)  =  IL(  IF) 

LW(K)  =  IW(IE) 

LH(K)  =  LH(K)  -  IH(IE) 

IF  (  LH(K)  .GE.  IH(IE)  )  GO  TO  128 
440   CONTINUE 
6  0   ML(K)  =  IL(I) 

MW(K)  =  IW(I) 

MH(K)  =  MH(K)  -  IH(I) 

N(MG)  =  N(MC)  -  1 

IF  (  MH(L)  .EQ.  0  )  L  =  L  +  1 

IF  (  M  .GE.  K  )  GO  TO  400 

M  =  M  +  1 

LA  =  LA  +  ML(K)  *  MW(K) 

KL(K)  =  ML(K) 

KW(K)  =  MW(K) 

IC  =  KL(K) 

ID  =  KW(K) 

MTX(IC,41)  =  MTX(ICt41)  +  1 

MTXf  ICtID)  =  MTX(ICtlD)  +  1 
400   CONTINUE 
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C    COMMENCE  LOADING 
6A0 


650 

810 
630 

840 

710 
660 

820 


780 

740 

700 
634 

667 
63  3 

670 


730 
760 
680 
770 


750 


NAC 
ISL 
IF  ( 
IF  ( 
ISL  = 
GO  TO 
LPA  = 
NAC  = 
NL  = 


=  0 

=  100 
MTX( ISLT41) 


ISL 
ISL 

640 
0 

NAC 
IAL 


EG.  1 
-  1 


+  1 


.GT.  0  )  GO  TO  650 
)  GO  TO  99 


KSL  =  ISL 
NW  =  IAW 
IF  (  NL  .EQ.  0  ) 
IF  (  NL  .LT.  KSL 
IF  (  MTX(KSL,41) 
IF  (  KSL  .EC.  1 
-  1 


KSL 
KSL  =  KSL 
GO  TO  630 
NL  =  NL  - 
=  KSL 

KSW 


GO  TO  800 
)  KSL  =  NL 
.GT.  0  )  GO  TO  840 
)  GO  TO  800 


KSL 


KSW 

IF  (  KSW  .GT.  NW 

IF  (  MTX(KSL,KSW 

IF  (  KSW  .EC.  1 

KSW  =  KSW  -  1 

GO  TO  660 

LPA  =  LPA  +  KSW 

MTXC KSL, KSW)  =  M 

MTX(KSL,  41)  =  M 

NW  --  NW  -  KSW 

IF  (  NW  .NE.  0  ) 

GO  TO  810 

IF  (  KSL  -EC.  1 

LSL  =  KSL  -  1 

JW  =  1 

JL  =  NW 

IF  (  MTX(LSL,41) 

IF  (  LSL  .EG.  1 

LSL  =  LSL  -  1 

JW  =  JW  +  1 

GO  TO  700 

DO  633  J  =  1,  NW 

IF  (  MTX(LSL,J) 

CONTINUE 

GO  TO  634 

IF  (  JL  .GE 

IDUM  =  JL 

JL  =  JW 

JW  =  IDUM 

JSL  =  ISL 

IF  (  JSL  .GT.  JL 

IF  (  MTX(JSL,41) 

IF  (  JSL  .EC  1 

JSL  =  JSL  -  1 

GO  TO  680 

LSW  =  LSL 


)  KSW  =  NW 
) .GT.  0  )  GO  TO  820 
)  GO  TO  780 


*  KSL 

TX(KSL,KSW) 
TX(KSL,  41) 

GO  TO  710 

)  GO  TO  800 


.GT.  0  )  GO  TO  667 
)  GO  TO  810 


.GT.  0  )  GO  TO  670 


JW  )  GO  TO  730 


)  JSL  =  JL 

.GT.  0  )  GO  TO  720 
)  GO  TO  750 
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790 
850 


870 

720 
880 

890 


860 

800 
900 

99 
57 


IF  (  LSW 
IF  (  MTX 
IF  (  LSW 
LSW  =  LS 
GO  TO  85 
IF  (  LSL 
LSL  =  LS 
GO  TO  74 
JSW  =  JS 
IF  (  JSW 


.GT. 
(LSL, 

.EQ. 
-  1 


NW  )  LSW  =  NW 
LSW).CT.  0  )  GO  TO  860 
1  )  GO  TO  870 


EC. 
-  1 


( 


MTX 

JSW 

=  JS 


IF 

IF  ( 
JSW 

GO  TO  88 
JL  =  JL 
MTX< JSL, 
MTX( JSLT 
LPA  =  LP 
IF  (  JL 
GO  TO  76 
NW  =  NW 
MTX( LSL, 
MTX(LSL, 
LPA  =  LP 
IF  (  NW 
GO  TO  79 
PLA  = 
AEF  = 
WRITE 
FORM  AT  (  « 
6  •  UNUSE 
7F4.2///) 
GO  TO  64 
AL  =  LA 
TVS  =  NV 
TAS  =  AL 
TTV  -  TV 
WRITE  (6 
FORMAT  < 

1  *  TOTAL 

2  ■  TOTAL 
STOP 
END 


LP 
PL 
(6 


W 
0 

L 
0 
L 

.GT. 
(  JSL, 

.EG. 
W  -  1 
0 

-  JSL 
41  ) 

JSW) 
A  +  J 
.EQ. 
0 

-  LSW 
LSW) 

41) 
A  +  L 
.EQ. 
0 
A 

A  /  ( 
,900) 

A/C 
D  LEN 

0 


1  )  GO  TO  810 


JW  )  JSW  =  JW 

JSW1.GT.  0  )  GO 

1  )  GO  TO  770 


=  MTX( JSL,  41) 
=  MTX(JSL,JSW) 
SW  *  JSL 
0  )  GO  TO  750 


=  MTX(LSL,LSW) 
=  MTX(LSL,  41) 
SW  *  LSL 
0  )  GO  TO  810 


TO  890 


IAL  *  I  AW  ) 

NAC,  LPA,  NL,  AEF 
NO.  •,  13,  5X,  •  AREA  LOADED  =  ',  15,  5X, 
GTH  =  •,  13,  5X,  •  AREA  EFFICIENCY  =  ', 


/  (  AL  *  IAH  ) 

/  (  NAC  *  IAL  *  I  AW  ) 

S  *  TAS 

,  57 >  TVS,  TAS,  TTV 
•  TOTAL  VERTICAL  EFFICIENCY    =  «,  F4.2,// 
AREA  EFFICIENCY        =  •,  F4.2// 
VOLUMETRIC  EFFICIENCY  =  »,  F4.2,///) 
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